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Abstract. We describe a search underway for periodic gravitational waves from 
the central compact object in the supernova remnant Cassiopeia A. The object is 
the youngest likely neutron star in the Galaxy. Its position is well known, but the 
object does not pulse in any electromagnetic radiation band and thus presents a 
challenge in searching the parameter space of frequency and frequency derivatives. 
We estimate that a fully coherent search can, with a reasonable amount of time 
on a computing cluster, achieve a sensitivity at which it is theoretically possible 
(though not likely) to observe a signal even with the initial LIGO noise spectrum. 
Cassiopeia A is only the second object after the Crab pulsar for which this is true. 
The search method described here can also obtain interesting results for similar 
objects with current LIGO sensitivity. 



PACS numbers: 04.80.Nn, 97.60.Bw, 97.60.Jd 
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1. Introduction 

The LIGO Scientific Collaboration (LSC) has so far published three types of searches 
for periodic gravitational waves (GWs): searches for known non-accreting pulsars 
[TJ [U [31 3], for the non-pulsing low-mass X-ray binary Sco X-l [5], and all-sky 
searches for as yet unknown neutron stars [SJ [3 E] . The first and last types of search 
are approaching the indirect upper limits on gravitational wave emission inferred from 
the observed spindowns (spin frequency derivatives) of pulsars and supernova-based 
estimates of the neutron star population of the galaxy [5] . 

Here we discuss the first of a fourth type of search for periodic gravitational waves: 
directed searches, which target likely neutron stars whose sky position is known to high 
accuracy, but whose spin frequencies and frequency evolution are not known at all. We 
describe such a search, which is currently underway, directed at the central compact 
object in the supernova remnant Cassiopeia A (Cas A). The data analysis challenge 
is to search a large parameter space of possible frequencies and frequency evolutions. 
We describe the object, estimate the computational costs of the search, and show 
that when the search of data from LIGO's recently completed S5 run is completed, it 
will beat the indirect limit on GW strain for Cas A. We also indicate how cost and 
sensitivity estimates can be extended to other directed searches. 

2. The central compact object in Cas A 

Cas A is a core-collapse supernova remnant, currently the youngest known in the 
Galaxy [10]. A central X-ray point source was discovered in first-light images taken 
by the Chandra X-Ray Observatory, indicating the presence of a central compact 
object (CCO). The nature of the CCO remains uncertain. No radio pulsations or 
7-ray emission have been observed, and there is no pulsar wind nebula observed in 
X-ray or radio; it is unlikely therefore that the CCO is an active pulsar [IT] . Proposed 
explanations include that it might be a young radio-quiet neutron star, or an accretion 
disk associated with a neutron star or black hole, or that it might be related to a type 
of slowly rotating neutron star known as an anomalous X-ray pulsar (AXP) or a 
soft 7-ray repeater (SGR) [UJH2]. Only in the first scenario could GW emission be 
detectable by LIGO. What makes Cas A an attractive target is its youth: the stars 
with the highest indirect limits (see next section) on gravitational radiation are young, 
and one could argue on theoretical grounds that any deformations left over from the 
violent birth of the star have had less time to be smoothed away by mechanisms 
such as viscoelastic creep. Young stars also spin more quickly than old ones. Of the 
seven confirmed CCOs, only two (possibly three) have measured spin periods [13] . 
The fastest is radiating gravitational waves at 20 Hz, just below the LIGO frequency 
band, but the other CCOs are also much older than Cas A. 

For the purpose of a directed search, we need to know the object's right ascension 
and declination. Chandra observations [11] have obtained these to sub-arcsecond 
accuracy [a = 23h 23m (27.945 ± 0.05)s, 5 = 58°48'(42.51 ± 0.4)"], which is sufficient 
for any GW observation. In order to define the range of search parameters and give 
an indirect limit on GW emission from the object, we also need the distance, age, and 
moment of inertia. The distance to Cas A has been estimated from the radial velocities 
of knots of ejected material to be 3.4+q ^ kpc [14] . Extrapolation of the proper motions 
of outer ejecta knots suggest a convergence date of 1681 ± 19, consistent with a possible 
observation by John Flamsteed in 1680 10J. Since computational costs are higher for 
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younger objects, we play it safe by taking 300 years (the approximate lower bound) as 
our fiducial age estimate. In what follows we use the canonical neutron star moment 
of inertia of 10 45 g cm 2 , although modern equations of state predict values higher for 
most neutron stars by a factor 2 or 3 [T5] . 



3. Indirect limits 



Indirect limits on the gravitational wave emission from rotating neutron stars are 
found by assuming that the gravitational wave luminosity is bounded by the time 
derivative of the total rotational kinetic energy: 

where e is the equatorial ellipticity, I zz the principal moment of inertia (assumed 
constant), and / the gravitational wave frequency (assumed to be twice the spin 
frequency) [5J [TB] . This condition is rearranged to give the "spindown" upper bounds 
on the ellipticity and the GW strain tensor amplitude Hq: 



5c 5 -/ 1 5GI zz -f 

e -V ^gT^' hQ -Dy—T' () 

The second limit is found from the first by substituting 
4n 2 G I zz ef 2 

h o = — 4 pp- ■ (3 

where D is the distance of the source [17] . 

For a directed search, the GW frequency / and its time derivative / are unknown, 
but the age is known. If we assume that the star is spinning down with f <x f n , and 
that it is currently spinning significantly more slowly than it was at birth, we can relate 
the frequency evolution to the characteristic age r and braking index n by [TH1 [THl [5] 

1 { f \ ff 

n = -rr. (4) 



n-l\-fj P 
If the spindown is dominated by GW from a constant mass quadrupole, then n = 5 
and t is the true age of the star. Substituting into the spindown limits (|2} gives 



T ftago < ■ ( 5 ) 



d8C " y l28ir 4 GI zz Tp ' * gL - D V 8c 3 r 

Using the numbers for Cas A from the previous section we get 



< 3.9 x 10- (^) 2 J( 10 7 Cm2 ) ( 3 °°; earS 



(7) 



Below we will consider searches over the range n = 2-7, including the possibility that 
n has changed since the supernova and thus a lifetime-averaged value is appropriate. 
Considering this, the uncertainty in D, and the fact that I zz may be triple our fiducial 
value (see discussion in (3|), these fiducial indirect upper limits are uncertain by about 
a factor of 2. Some theories of quark matter allow for ellipticities in the range indicated, 
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though normal neutron star models do not [201 US HI] ■ An internal magnetic field of 
order 10 16 G could also produce such ellipticities [33J 12H HB] , although it is not 
clear if such a field is stable, and if the external field is this strong then the star by 
now has spun down out of the LIGO frequency band. The age-based indirect limits 
serve, like the spindown limits, as indicators of which objects are interesting, but since 
they are based on less information they are not as solid as the spindown limits. It is 
not known if Cas A spins in the LIGO band (period < 50 ms), and indeed only 10% 
of known pulsars do so [27]. Thus a search such as we describe could detect an object 
on the speculative end of the range of theoretical predictions. 

4. Search method 

The LSC uses both fully coherent [2 H USUI] and semi - coherent [3 El 09 US methods to 
search for periodic gravitational waves. Semi-coherent methods are computationally 
cheaper than coherent methods, but coherent methods can achieve greater sensitivity 
if the cost is feasible. 

For a young neutron star such as Cas A the integration time needed is short 
enough (see next section) for us to pursue enhanced sensitivity without undue 
computational cost. We therefore use the fully coherent JF-statistic search [5], as 
implemented by the ComputeFStatistic_v2 routine in the LSC Algorithm Library 
|28j . This routine computes optimal filters for the gravitational wave signal, including 
modulation by the detector beam patterns, in multiple interferometers which are 
treated as a coherent network [T71 [55] . This search uses data from the 4km LIGO 
interferometers at Hanford, WA, and Livingston, LA. 

The computation is conducted in the frequency domain using short Fourier 
transforms (SFTs) of segments of strain data, typically of 30 minutes duration so 
that the GW frequency will remain in one frequency bin over the length of the SFT 
[5]. The SFTs are vetoed by a suite of data quality flags to remove poorer quality 
data. For windows of up to 15 days during the first year of the S5 run the duty cycle - 
the ratio of post-veto SFT live time to total time span, averaged over interferometers 
- can somewhat exceed 70%. 

A search for a young neutron star such as Cas A, which is younger than objects 
considered in previous LIGO multi-template searches, must cover a greater spindown 
parameter space including a second frequency derivative (see next section). This has 
required the extension of existing LSC software to efficiently cover a three-dimensional 
space using the parameter space metric. The points are distributed on a body-centered 
cubic (bec or A^) lattice, which is known to be the optimal lattice covering in three 
dimensions [30] . 

In the event no plausible signal is found, we will set upper limits by methods 
similar to the frequentist analyses in [TJ [S]. These are based on Monte Carlo 
simulations searching the data for a multitude of software-injected signals with 
a distribution of amplitudes, inclination angles, and polarization angles in each 
frequency bin. We will also test on a smaller set of simulated signals which were 
hardware injected into the S5 data. 
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5. Estimated cost and sensitivity 

The sensitivity of a search for periodic signals can be put in terms of the 95% confidence 
limit on GW strain tensor amplitude, which takes the form 

hf % = Q^S h {f)/T Aat . (8) 

Here Sh is the strain noise power spectral density, Td a t is the data live time, and Q is a 
statistical threshold factor which depends on the parameter space and other details of 
the data analysis pipeline. For a coherent multi-interferometer search, the limits add 
in inverse quadrature. Monte Carlo simulations searching for injected signals from 
Cas A, as well as the results of the similar multi-template Crab search [4], indicate 
that is in the mid-30s for a directed search, and thus we use 35 in our estimates 
below. Because is determined by the tail of a Gaussian distribution, it is very 
weakly dependent on the volume of parameter space searched. However the data live 
time Tdat is computationally limited and thus does depend on the parameter space. 

The parameter space range is chosen as follows. The frequency band is chosen to 
be 100-300 Hz, which surrounds the band where the LIGO interferometers are most 
sensitive. As we shall see below, this is roughly the band over which a directed search 
can beat the indirect limit on ho with reasonable computational cost. The frequency 
derivative ranges are chosen based on considering braking indices n in the range 2-7. 
This range covers all known pulsars, except the Vela pulsar which is visibly interacting 
with its wind nebula (nonexistent for Cas A). It also includes the values for radiation 
dominated by a static dipole or quadrupole (n = 3 or 5) as well as a saturated r- 
modc (n = 7) [31] . Thus the range of each frequency derivative depends on the lower 
derivatives, and we have 

100 Hz < / < 300 Hz, L<-f<L t 2 Jl<f < III. (9) 
6r t f f 

Note that the range of / by definition is related to the present-day braking index, 
while the range of / corresponds to an average braking index over the lifetime of the 
star. Thus we allow for the braking index varying over time between the indicated 
limits. 

There remains the problem of efficiently tiling, or choosing specific points in 
parameter space for which to compute the ^-statistic. It is straightforward to apply 
the method of [32] to find the parameter space metric [33J 

4^1n +2 (j + l)(fc+l) 

l3k (j' + 2)!(fc + 2)!Cj + * + 3)' [W > 

where the components are with respect to the fcth derivative of the GW frequency 
at the beginning of the observation, T span is the total duration of data (including 
dropouts), and the indices j, k take integer values between and the highest derivative 
considered (2 for Cas A). This metric, which is the Fisher information matrix with a 
phase constant projected out, is used to set up an efficient tiling which takes advantage 
of the covariances between parameters. The number of points needed for an optimal 
(bec or Al) tiling is given by [34] 

N p ~ 0.19/i~ 3/2 A /det^ (11) 

T 

where fi is the mismatch and we have performed the integral in equation (24) of [34] 
using the ranges ([9]) and discarding the lower bound on frequency, which is only a few 
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Figure 1. Estimated sensitivity of an S5 search compared to the indirect limit 
on GW emission for Cas A. Search parameters are the fiducial ones described in 
the text. 



percent correction. We determine the highest frequency derivative needed by finding 
k such that 7fcfcAj: > fi, where is the range of the fcth frequency derivative and 
we take /j, to be 20% (typical for periodic signal searches). In our case / is required 
for T span greater than about a week; as shown below, this applies for any search 
competitive with the indirect limit. 

Since equation (jlip is obtained by dividing the proper volume of the parameter 
space by the proper volume per template, we expect it to underestimate N p of a 
practical implementation due to the need to cover the edges of the parameter space. 
Because the extent of our parameter space in / is often comparable to or less than the 
unit cell length of a single template, we expect that an ideal lattice covering would 
require several times the ideal number of templates in (jlip . Technical limitations of a 
speedy — and therefore simple — template bank generation algorithm also require us to 
lay extra templates to guarantee that the edges of the parameter space are completely 
covered. We have found from Monte Carlo simulations that the combination of these 
effects can cause (jlip to underestimate N p by up to an order of magnitude. Even 
in this worst case, without any improvement of existing template bank algorithms, 
the computational cost is still feasible since our fiducial estimate below is for a small 
number of computing nodes. The size of the template bank should not significantly 
affect the upper limits, which are very weakly dependent on the number of templates 
and thus on the number of statistical trials. 

Finally we estimate the computational cost and sensitivity of a directed search. 
Preliminary runs on nodes of the APAC cluster [35] find a timing of about 6 x 10~ 7 s 
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per template per SFT. Assuming 30-minute SFTs and two interferometers with 70% 
duty cycle, the computing time for the search (exclusive of Monte Carlo simulations 
to compute upper limits) is 

/max V ^300 years ^ 3 ( T span \ 7 f 200 



20d " ys [mk) [=^) [w^;) \j£s)- ( 12 » 

For these fiducial parameters and two interferometers with the initial LIGO design 
noise spectrum (33] and 70% duty cycle, the sensitivity curve flSJ) is plotted in figure H 
The minimum of the curve (smallest detectable ho) is 

where the latter scalings allow T span to vary at fixed computational cost and are 
useful for evaluating searches for other objects. Combining the previous two equations 
indicates that the sensitivity only improves as the 14th root of the computational cost, 
and thus there is not much point in integrating for significantly longer without an 
improved semi-coherent analysis method. 

Thus we see that this search for Cas A, when completed on S5 data, will beat the 
fiducial indirect limit on GW emission from about 100 to 300 Hz. This will double 
the number of objects (after the Crab pulsar) for which initial LIGO has beaten an 
indirect limit. Similar searches can be made for other suspected young neutron stars. 
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